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Abstract 

Shear fields due to weak gravitational lensing have characteristic coherent patterns. We describe 
the topological defects in shear fields in terms of the curvature of the surface described by the lensing 
potential. A simple interpretation of the characteristic defects is given in terms of the the umbilical 
points of the potential surface produced by ellipsoidal halos. We show simulated lensing shear maps 
and point out the typical defect configurations. Finally, we show how statistical properties such 
as the abundance of defects can be expressed in terms of the correlation function of the lensing 
potential. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 



The central tenet of Einstein's General Relativity is that massive bodies curve space- 
time. As a result, light from distant galaxies is deflected by mass distributions encountered 
along the line of sight. The images of distant galaxies, that act as sources, are magnified 
and sheared - this effect is known as gravitational lensing. The most striking manifestation 
of gravitational lensing effects, know as strong lensing, consists in the formation of multiple 
images of a single background galaxy. In the special case in which a source, a very large 
mass and the observer happen to be approximately aligned an Einstein ring can be observed 
and from its diameter the lensing mass can be inferred 

Along typical lines of sight the lensing effect is weak, leading to percent level magnifi- 
cations and shears. However, such a small signal can still be detected from a statistical 
analysis that relies on the coherence of the shear field over the sky. This allows us to infer 
how the dark matter is concentrated around galaxies and galaxy clusters, as well as providing 
a testing ground for dark energy and modified gravity theories {2, sl. 

In this paper we explore the connection between the theory of topological defects and 
the spatial patterns of shear fields due to weak gravitational lensing. The starting point 
of our approach rests on an analogy between gravitational lensing shear fields, as a probe 
of structure formation on cosmological scales, and the anisotropic optical or mechanical 
response of materials, as a probe of their inhomogeneous structure on microscopic scales. As 
an illustration, the topological defects in the local shear field of an elastic medium reflect the 
external deformations applied to the solid [4]. Similarly for thin liquid crystal films confined 
on a curved substrate, the density of topological defects depends on the inhomogeneous 
curvature of the underlying surface [s],!^, 3]- In a somewhat different context, the distribution 
of optical singularities can be used to shed light on the statistics of randomly polarized light 
fields ^] . Here we suggest that topological defects in the cosmic shear can be used as a probe 
of the gravitational potential generated by the lensing mass fluctuations on large scales. 

The outline of this work is as follows. In §2, the basic formalism of weak gravitational 
lensing and the criteria to identify topological defects in shear fields are presented. Our 
geometric approach is presented in §3 where the defects in the shear field are related to 
the properties of an imaginary surface whose lines of constant height are defined by the 
contour lines of the underlying gravitational potential responsible for the lensing. This 
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mapping is applied to describe the characteristic behavior of eUipsoidal structures in the 
mass distribution and to more reahstic shear lensing maps obtained using N-body computer 
simulations. In §4 we turn to the case of a random mass distribution whose gravitational 
potential is a Gaussian variable that describes the stochastic geometry of a surface. This 
assumption (valid on large angular scales) allows us to express the defect density in terms of 
the two point correlation function of the gravitational potential. The latter is estimated for 
the standard Cold Dark Matter model for large-scale structure in the universe and compared 
to the results of simulations. Corrections for weakly non-Gaussian fields can be calculated 
using perturbation theory as discussed in §5 where ideas for further work are briefly sketched. 



II. TOPOLOGICAL DEFECTS IN LENSING SHEAR FIELDS 

In this section, we review the basic formalism of gravitational lensing that relates the 
deformation of the shapes of background galaxy images to the projected mass density re- 
sponsible for lensing. 



A. Basics of Gravitational Lensing 

Consider a source whose true angular position on the sky makes an angle j3 with an 
arbitrary optic axis. As a result of the lensing, an observer sees the light ray as coming 
from an image at an angle 9 (from the optic axis) that differs from (3 by the (reduced) 
deflection angle a. These angular displacements can also be viewed as two dimensional 
vectors {a, (3,6} on a locally fiat sky (see [Oj for a review). These vectors are related by: 
^=6 -a. 

In this notation, the image distortion caused by gravitational lensing is locally described 
by Aij, the Jacobian matrix of the transformation between (3 and 9 which reads 

- 09, ~ 89, ~ d9,d9j 

where we used the relation a = Vg- between the deflection angle and the gradient of 
the projected (two-dimensional) Newtonian potential ip{9). The latter satisfies the Poisson 
equation 

li^U + i^22) ^ 1^1^ = m (2) 



tion field, see Ref . m for precise definitions. Equation ([2]) fixes the trace of the deformation 
matrix Aij. The other two independent components of Aij can be rewritten in terms of a 
shear tensor 7 = {71,72} and the angle a as 

71 W = ^(^ii-V'22)=7Wcos[2a(^"*) 

72(^"') = ^-12 = ^-21 = 7(^1 sin [2a{e)j (3) 

1 /2 

where 7 = (7^ + 7I) > and a is the orientation angle of 7. The physical interpretation 
of the shear tensor is straightforward once the components of the Hessian ipij are expressed 
in terms of 71, 72 and n with the aid of Equations ([2113]) and then substituted into Eq. (Q. 
The result reads 

(1 \ / cos 2a sin 2a \ 
-7 (4) 
1/ \ sin 2a — cos 2a I 

The eigenvalues of this matrix are {1 — k — 7,! — K + 7} with eigenvectors vi = [cos a, sin a] 
and V2 = [— sin a, cos a], respectively. The direction of shear is vi: when the deformation 
tensor {A~^)ij acts on a test circular source, it magnifies the image isotropically, by the first 
term in Eq. and it deforms it into an ellipse with major axis of magnitude (1 — k — 7)"^ 
oriented in the a direction, by the second term (see Fig. [T]). 

In the sub-section below we will consider the pattern of the shear on the sky induced 
by the spatial variation of ip{6). Over patches of order a degree on a side, characteristic 
tangential patterns around lensing galaxy clusters and galaxy groups are observable. 

Thus measuring the anisotropy introduced in the light distribution from distant galax- 
ies allows us to infer information about mass inhomogeneities along the line of sight that 
generated the lensing shear. 



B. Topological defects 

The local orientation of the shear can be visualized by means of a line field, that is to 
say a vector n with both ends identified n = — n. This is analogous to the description of 
interacting liquid crystal molecules in condensed matter physics, where the director field n 
denotes the local orientation of the molecules. In that context, it is often useful to develop 
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FIG. 1: The light from a distant source of circular shape (drawn as a dashed line) is deformed into 
an ellipse (continuous line). The local direction of the major axis of the ellipse can be represented 
by a double headed vector that forms an angle a with the x axis of a cartesian coordinate system. 
The positions in a small patch of the sky are measured by the angles {^1,^2} within the flat sky 
approximation . 

a coarse grained view of the many-body system by concentrating on global configurations of 



the director field called disclinations Disclinations are topological defects like vortices 
in fluid mechanics that can be classified according to the index of the streamlines around 
them. 

For example , the shear field angle a{9i, 62) has a topological defect of index m if 

Va{ei,e2) ■ dl = 2Tim (5) 



around any contour that surrounds the location of the defect 6^. Upon using Stokes theorem, 
Eq. ([5]) can be cast in the differential form 

V X Va(0i, 62) = 27Tm 61(9- dd) (6) 

where 5j~,{9 — 9d) denotes the Dirac delta function at the location of the defect. As a result, 
topological defects are analogous to localized regions of quantized magnetic flux with charge 
proportional to the winding number m, and Va plays the role of the electromagnetic gauge 
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FIG. 2: (Color Online) Basic classification of topological defects labeled by their respective index. 
The index m is found by evaluating a line integral of Va around any contour surrounding the 
singularity at the core of the topological defect and dividing by 2tt, see Equations ([T]) and ([2]). 

field. The (shear) field lines corresponding to topological defects of different indexes are 
reproduced schematically in Fig. [21 The symmetry of the shear field requires that the defect 
index m must be a half-integer. The negatively charged defects correspond to the case in 
which the shear field changes in the clockwise direction while travelling counter-clockwise 
around the contour encircling the defect. 

Global rotations by a constant angle at every point in the fiat sky, a{6) ^ a{6) + c, do 
not change the index of a topological defect since the definition in Eq. (jl]) depends only 
on the gradient of a. For example, the vortex-like defect in Fig. [2] has the same m = 1 
index as a sink which is obtained from the vortex by a c = 7r/2 rotation. This invariance 
distinguishes the topological classification adopted in this work from the more common E- 
B modes classification adopted in previous studies of the CMB polarization field and of 
gravitational lensing. For example a 45 degree rotation of each shear around an m = 1 
defect does not change its topological charge, but it changes an E mode to a B mode. In 
this study, we concentrate on global topological features of the shear field whose positions 
on the sky are related (non-locally) to the large scale mass fluctuations responsible for the 
lensing. 
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III. THE INDUCED GEOMETRY OF THE LENSING SHEAR 



In this section, we derive the connection that exists between the topological defects 
in the shear field and the projected gravitational potential ip{6) generated by the mass 
fiuctuations. It is convenient to view ip{9i,92) as the height function of a non-intersecting 
two dimensional surface. In the weak lensing regime, the derivatives of the potential ip are 
small; the differential geometry of the surface is then approximated by the properties of the 
Hessian ipij [12|. We show that the shear direction corresponds to the principal direction of 
maximal curvature on the surface, provided that 7 ^ 1. 

A. Defects as umbilical points 

On such a two dimensional surface, we can step away from any given point in an infinite 
number of directions so that for each direction a curvature is defined. As the direction is 
smoothly varied, two perpendicular directions of principal curvatures can be found for which 
the curvature is maximal and minimal. The two principal curvatures, henceforth denoted 
by {/ti, fi;2}, are the eigenvalues of the Hessian matrix ipij-. 

{/ti, K2} = i (^tpu + ^22 ± ■\/4^?2 + (^11 - ^22)'^ = K ± 7 (7) 

The eigenvectors give the local directions of the principal axes. Since the magnitude of 
the shear field 7 is positive, ki = k + j corresponds to the principal direction of maximal 
curvature. The Hessian matrix ipij determines the deformation matrix Aij according to Eq. 
([1]). Upon comparing Eq. ([7]) with Eq. (jl]) and the discussion following it, we can conclude 
that the direction of maximal curvature is the shear direction along which a reference circular 
source is stretched most by the lensing potential (eg. it points along the major axis of the 
ellipse in Fig. [1]). 

The two principal curvatures are equal when 4-?/'^2 + (V'li ~ "^22)^ is equal to zero. This 
defines the umbilical points of the surface: points where the principal directions are undefined 
and the surface is locally spherical or fiat fl3|. 

We can now connect the geometry of the surface which represents the gravitational po- 
tential and the topological defects in the lensing shear field. At points where the shear field 
7 has a topological defect, both of its components, see Eq. ([3]) must vanish, because its 
local direction is undetermined. This ensures that the two principal curvatures are equal. 
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As a result, topological defects of the shear field can be identified as the umbilical points of 
the induced surface. Indeed the umbilical points can be classified according to the index of 
the principal direction vector field, which is either +i or — | at an umbilical point. As we 
shall see this corresponds to topological defects of index ±| in the shear field. 

The local angle a{6i, 62) by which the coordinate axis needs to be rotated to overlap with 
the principal axes (given by the eigenvectors of the Hessian of the height function ip) is 

tan 2a = , (8) 
- V22 

A comparison of Eq. ([8]) with the definition of the shear components 71 and 72 in Eq. ([3]) 
shows that a = a, the angle that specifies the orientation of the shear 7. Thus the principal 
direction of maximal curvature on the surface tracks the shear field. 

In the next sections we will show how this mathematical mapping leads to connections 
between the defect locations and the distribution of mass that generates the gravitational 
potential. In addition the identification between topological defects and umbilical points of 
a random surface will allow us to relate the density of defects to statistical quantities such 
as the two-point correlation function of the projected gravitational field. 

In practice, one would attempt to identify defects either in direct experimental measure- 
ments of shear using galaxy images, or ray tracing computer simulations that provide the 
shear field on a discrete grid. In both cases, it is necessary to consider how to interpret the 
notion of a defect (introduced via Eq. ([5]) in the continuum limit) on a square grid. We have 
carried out such a measurement using the net change in the angle a around closed loops at 

n 

each pixel vertex, using techniques similar to those of p^]. Further details of the numerical 
aspects will be presented elsewhere. 



B. The ellipticity of haloes 

To illustrate the formal ideas presented in the last section consider the simple case of an 
elliptical potential field ip{x,y) generated, for example, by an isolated distribution of mass 
which is not axisymmetric, namely an elliptical halo. (Note, however, that even if the mass 
distribution is locally axisymmetric the corresponding potential can be perturbed into an 
elliptical shape by the tidal field of nearby objects, see next section [15].). 

For simplicity, consider two simple models of elliptical haloes whose gravitational poten- 
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FIG. 3: Panel (a) and (b) show contour plots of the gravitational potentials for the elliptical halos 
in Eq. Q and Eq. (jlOp respectively (dark blue is most negative). The parameters were set to 
e = 0.3 and c = 0.2. The two black dots indicate the umbilical points that correspond to the two 
+1/2 topological defects in the shear field. The wiskers represent the local direction of the cosmic 
shear which is oriented along the principal direction with maximum curvature. 



tials are given in dimensionless variables by 

ip{9i,92) = -exp 

and 



01 



1 + e 



+ 



61 



(9) 



(10) 



Figures [3]^a) and MJo) show the contours plots of the two potentials in Eqs. and ffTU]) re- 
spectively, overlayed onto the corresponding (normalized) shear fields represented by wiskers. 
The umbilical points are labelled by black dots. Comparison of the shear field around each 
defect in this figure with the defects shown in Fig. [2] allows a clear identification of the 
two umbilical points with two +1/2 topological defects. The line joining the two defects is 
perpendicular to the major axis of the ellipsoidal contours of the potential (and of the mass). 
This follows from the fact that the shear field is oriented along the principal direction of 
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maximal curvature. The defects separation, s, measured in the same units of the angular 
distances {01,62} adopted for the potential in Eq. Qj, reads 



1 + e 

while the separation for the potential in Eq. ffTOj) reads 



(11) 



In the limit of zero ellipticity, e = and the separation between the two defects, s ~ ^/e, 
goes to zero. The two +1/2 defects coalesce into a single vortex of index +1 which is 
expected from an isolated axisymmetric mass distribution [9|. Qualitatively similar behavior 
is expected for ellipsoidal haloes with different potential profiles. We have also checked this 
for a non-singular logarithmic potential. 

Hence, the separation of the two nearby defects is a measure of the ellipticity of the 
haloes. This is potentially useful since large-scale structure studies have shown that the mass 
distribution in the universe can be well approximated (especially for statistical purposes) as 
a network of elliptical halos of varying mass and concentration . Thus our approach 

facilitates the interpretation of observed shear fields, as we discuss further below. 

We note that the analysis presented here applies to non-singular mass distributions which 
generate a gravitational potential that is both continuous and finite everywhere, so that 
the induced surface is smooth. In the case of either a point mass or an isothermal mass 
distribution, the lensing potential diverges at the origin and would require further refined 
techniques. We expect that realistic mass profiles and the ellipsoidal nature of halos avoids 
such a divergence. 



C. Complex defects patterns in realistic shear maps 

The geometric analysis developed in this work can be applied to shear maps generated by 
more realistic mass distributions than the isolated halo considered in the previous section. 
For instance, tidal effects between nearby mass distributions can generate complex defect 
patterns that capture the skeleton of the contour lines of the corresponding gravitational 
potential. In order to demonstrate this point, we have generated shear lensing maps by 
means of N-body simulations combined with a ray tracing algorithm 
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Figure H] shows a patch of the observed sky with the simulated (normahzed) shear field 
plotted as whiskers at each pixel with the overlaid color map of the projected mass distribu- 
tion (contours of k). The continuous line interpolates the whiskers to facilitate comparison 
with the defect patterns reproduced in Fig. [21 The red and black dots indicate negative 
or positive defects whose index has magnitude 1/2 (top row in Fig. [2]). They have been 
identified by measuring the net change, Aa, accumulated by the shear angle in going around 



each pixel, see Ref. 



The top left part of the plot shows a prominent mass over-density of irregular shape 
which generates at large scales shear field lines topologically equivalent to a vortex. Upon 
magnification, two +1/2 defects can be resolved. The line along which the two defects lie 
is approximately oriented perpendicularly to the major axis of the nearly elliptical mass 
distribution, as indicated by our geometric analysis in the previous section. 

In addition, note the presence of the isolated —1/2 defect (red dot) generated by tidal 
effects at the center of three mass over densities. The local projected gravitational potential 
resembles a monkey saddle when plotted as a height function with an umbilical point (the - 
1/2 defect) at its center. Note, in addition, that the dipole of ±1/2 defects in close proximity 
(on the right side of the picture) has no significant effect on the field lines at large scales 
since the two defects cancel each other. 



IV. STOCHASTIC GEOMETRY AND MASS FLUCTUATIONS 

In this section, we consider the case of a random mass distribution and present a field 
theory from which the density of topological defects can be read off from the statistics of 
the stationary random function ilj{x,y). The gravitational potential is first assumed to be a 
Gaussian random variable (deviations from Gaussianity can be calculated using perturbation 
theory, this is beyond the scope of this paper). Light coming from very distant sources 
undergoes deflections from a number of intervening lensing planes. Assuming that the mass 
distribution in different redshift slices is uncorrelated, the projected gravitational potential 
resulting from a large number of such lensing planes is well described by a Gaussian random 
variable according to the central limit theorem. In addition, with a large smoothing scale, the 
lensing mass distribution itself is well approximated by a Gaussian random field. In practice, 
there are not sufficiently many uncorrelated lens planes and the Gaussian assumption is valid 
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FIG. 4: Shear lensing map generated from ray tracing simulations on a patch of the sky with 
angular size around 4 arcminutes. The wiskers represent the local direction of the shear field 
interpolated by the continuous line. The colour plot is the projected mass distribution and the 
defects of index plus or minus 1/2 are indicated by black and red dots respectively. 

only on angular scales sufficiently larger than 10 arcminutes (it is a better approximation 
for sources at higher redshift). 



A. Defect density 



In order to count A^, the number of topological defects (where the magintude of 7 = 0) 
in a patch of the sky, we need to evaluate the surface integral 

d{iJn - ^"22, ^12) 



N 




dOi (192 (^d(^11 - ^22) ^Z)(^12) 



det 



(13) 



5(^1,^2) 

where 5d indicates a delta function and the appropriate Jacobian determinant for the change 



of variable has been inserted 



20[ |. This integral can be performed explicitly once the grav- 



itational potential ip{0i,62) and its derivatives are specified. Note that the third spatial 
derivatives are also needed to evaluate the determinant in Eq. (fT3l) . 

In the case of Gaussian fluctuations of the lensing potential, the average defect density 
can be obtained by performing a functional integral that averages over the unknown fields 
with a probability density which is simply given by the exponential of a quadratic function 
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of if) and its derivatives. Following the seminal work of Berry and Hannay 20], we can write 
down explicit formulas for the density and other statistical properties of the detects using 
standard field theoretic manipulations. The statistics of ip{9) are completely determined by 
specifying either the autocorrelation function C{9) or its power spectrum P{k) defined by 



C{e) = {ij{9o)ij{9o + 9))= / £kPik)e 



Jk-9 



(14) 

The number density of defects, d, can then be related to the autocorrelation function of the 
lensing potential 



d 



3C, 



(6) 



IOttC, 



(4) 



(15) 



where the derivative C^^^ indicates the sixth spatial derivative of the correlation function 
C{9) evaluated at ^ = ^bl]. 

This formula can also be cast in terms of the moments of the power spectrum P{k) from 
Eq. (iMl) . The result reads 

(16) 



d 



AttMa 



where the n moment of the power spectrum is defined as 



27T I dkk"+^ 

'o 



P{k) 



(17) 



The two correlation functions of topological charges with equal or opposite sign, 

{5'++(6') = g (9), g^^{9)}, can be readily obtained from Eq. f|T3|) in terms of 9, the angular 

distance in the projected sky between two defects of equal or opposite index. The resulting 
mathematical expressions are too complicated to list here, see 2m\ and references therein 
for more details. The topological defects correlation functions depend on the full functional 
form of C{9) and its derivatives, not only its asymptotic value for 9^0. As a result, they 
are a more sensitive probe of the underlying cosmological processes than the defect density 
in Eq. (fTSD. 

The mathematical techniques adopted to study the topology of the shear field can be 
applied successfully to the study of the cosmic microwave polarization field 1^, l22|, |23| . The 
components of the lensing shear field are always correlated because in the standard weak 
lensing approximation they are derived, via Eq. ([3]), from the Hessian of the gravitational 
potential which is the fundamental physical field that controls the statistics. (Departures 
from this behavior are indicative of systematic errors in the data or physical effects distinct 
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from lensing. ) The CMB polarization components have a B-mode contribution from pri- 
mordial tensor mode fluctuations and from lensing along the line-o f- sig ht, so one generally 
requires both a scalar and pseudoscalar potential to describe them [24 1. 



B. Application to the cold dark matter model 

The relations in Equations ( |T5lfT6i) express the number density of defects d in terms of 
the correlation function or power spectrum of the projected gravitational potential. We can 
evaluate these for the current cosmological best fit A-CDM model, with the caveat that 
the Gaussian approximation breaks down on scales below about 10 arcminutes (or angu- 
lar wavenumbers £ above about 1000), the precise value of the angle requires careful tests 
and depends on the source redshift. The breakdown occurs due to nonlinear gravitational 
dynamics which couples the initially random Fourier modes that comprise the perturbed 
potential [l^. While nonlinear effects are evident in measures like the lensing power spec- 
trum on larger scales, the Gaussian approximation may still be valid to study topological 
features which can be altered only by orbit crossings or mergers. 

The potential is obtained from the projected density field via the Poisson equation as 
discussed in §11. The density field is characterized by its power spectrum, which can be 
approximated analytically by 

P{k) ^ AkT\k) (18) 

where the transfer function T{k) can be approximated analytically in terms of g = k/Q^h'^ 
(Bardeen et al 1984) 

T{q) - ^^^^^"^3^^^^^^ [1 + 3.89g + (16.1g)^ + (5.46^)=^ + (6.71g)^] (19) 
This power spectrum, evolved in time using linear perturbation theory, and projected 



appropriately along the line of sight gives C{6), the two-point correlation function of ip 18|. 
Upon evaluating Eq. ( fTSl) above, we find (for comparison with simulations, evaluating the 
expression at an arcminute), the predicted number density of defects is ~ 10'' per steradian. 
The two main approximations made are the assumption of Gaussianity/linear evolution and 
the assumption that the asymptotic behavior as 6 approaches is recovered at scales of 
order an arcminute. 
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We compare the analytical prediction made above with the measured number of ±1/2 
defects in ten independent ~ 6 square degree patches obtained from simulated shear maps. 
We estimate the measured density is 2 x 10^. (The error estimated from the standard 
deviation is less than 5% .). We may expect that non-Gaussianity, due to the mergers of 
several small halos into larger ones, tends to reduce the number of defects. Hence, a difference 
within a factor of order unity is perhaps not surprising. A more careful comparison will be 
attempted by making perturbative corrections to the Gaussian formula in future work. 

V. CONCLUSION 

In this paper we have studied topological defects in the shear field generated from weak 
gravitational lensing. Our geometric approach rests on the observation that topological 
defects correspond to umbilical point of an imaginary surface whose height function is given 
by the projected gravitational potential. This allows us, for example, to relate the ellipticity 
of a gravitational halo to the distance between two defects of index +1/2. We describe the 
overall pattern of the shear field in terms of the defects generated by clustered halos, as 
shown in Figure 4. Moreover the density of the defects can yield information on the two 
point correlation function of the gravitational potential if the fluctuations are assumed to be 
Gaussian. In this case, the statistical properties of the shear field can be readily calculated 
from a simple field theory. 

Topological defects in the shear field provide a different view than direct measurements 
of shear correlations, which is the standard approach in lensing cosmology. Whether the 
study of defects can provide cosmological information is an open question. We have not 
studied the robustness of defect identification or their statistical properties in the presence 
of measurement noise. This should be straightforward to carry out, at least with Gaussian 
noise added to simulated shear fields. 

The analytical results in Section IV rely on the statistics of a Gaussian random field. On 
small scales, the lensing shear has distinct non-Gaussian features, so the Gaussian description 
is valid only on large enough scales or alternatively from sources at very high redshift. It 
may however be a suitable starting point for perturbative calculations of the effects of weak 
non-Gaussianities associated with the onset of the non-linear gravitational dynamics regime. 
In general, measurements from simulations can be used to test analytical results, to compare 
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with data and to find signatures of non-Gaussian features in the defect statistics. 
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